* use Helmert transform for OSBG36 <-> WGS84.
* adjust unicsv_grids reference for Helmert xform.
* update copyright.
* use constexpr to select helmert inversion method.
* tweak helmert.
#include "formspec.h" // for FormatSpecificDataList
#include "garmin_fs.h" // for garmin_fs_t
#include "garmin_tables.h" // for gt_display_modes_e, gt_find_desc_from_icon_number, gt_find_icon_number_from_desc, gt_get_mps_grid_longname, gt_lookup_datum_index, gt_lookup_grid_type, GDB, gt_get_icao_cc, gt_get_icao_country, gt_get_mps_datum_name, gt_waypt_class_names, GT_DISPLAY_MODE...
-#include "jeeps/gpsmath.h" // for GPS_Math_Known_Datum_To_UTM_EN, GPS_Math_WGS84_To_Known_Datum_M, GPS_Math_WGS84_To_Swiss_EN, GPS_Math_WGS84_To_UKOSMap_M
+#include "jeeps/gpsmath.h" // for GPS_Math_Known_Datum_To_UTM_EN, GPS_Math_WGS84_To_Known_Datum_M, GPS_Math_WGS84_To_Swiss_EN, GPS_Math_WGS84_To_UKOSMap_H
#include "src/core/datetime.h" // for DateTime
#include "src/core/logging.h" // for FatalMsg
#include "src/core/textstream.h" // for TextStream
case grid_bng:
- valid = GPS_Math_WGS84_To_UKOSMap_M(wpt->latitude, wpt->longitude, &east, &north, map);
+ valid = GPS_Math_WGS84_To_UKOSMap_H(wpt->latitude, wpt->longitude, &east, &north, map);
if (valid) {
*fout << QString::asprintf("%s %5.0f %5.0f\t", map, east, north);
}
** @author Copyright (C) 1999 Alan Bleasby
** @version 1.0
** @modified Dec 28 1999 Alan Bleasby. First version
+** @modified Copyright (C) 2024 Robert Lipe
** @@
**
** This library is free software; you can redistribute it and/or
********************************************************************/
#include "jeeps/gpsmath.h"
-#include <cassert> // for assert
#include <cmath> // for sin, tan, cos, pow, log, sqrt, asin, atan, exp, fabs, round
#include <cstdint> // for int32_t
#include <cstdlib> // for abs
#include "defs.h" // for case_ignore_strcmp, fatal, CSTR
#include "jeeps/gpsdatum.h" // for GPS_ODatum, GPS_OEllipse, GPS_Datums, GPS_Ellipses, UKNG, GPS_SDatum_Alias, GPS_SDatum, GPS_DatumAliases, GPS_PDatum, GPS_PDatum_Alias
+static constexpr bool use_exact_helmert_inverse = false;
static int32_t GPS_Math_LatLon_To_UTM_Param(double lat, double lon, int32_t* zone,
char* zc, double* Mc, double* E0,
}
+/* @func GPS_Math_Helmert *******************************************
+**
+** Transform one datum to another
+**
+** @param [r] Sx [double] source X
+** @param [r] Sy [double] source Y
+** @param [r] Sz [double] source Z
+** @param [w] Dx [double*] converted X
+** @param [w] Dy [double*] converted Y
+** @param [w] Dz [double*] converted Z
+** @param [r] tX [double] translation along x axis (meters)
+** @param [r] tY [double] translation along y axis (meters)
+** @param [r] tZ [double] translation along z axis (meters)
+** @param [r] sppm [double] scale factor - 1 (ppm)
+** @param [r] rXs [double] rotation about x axis (seconds)
+** @param [r] rYs [double] rotation about y axis (seconds)
+** @param [r] rZs [double] rotation about z axis (seconds)
+**
+** @return [void]
+************************************************************************/
+void GPS_Math_Helmert(double Sx, double Sy, double Sz,
+ double* Dx, double* Dy, double* Dz,
+ double tX, double tY, double tZ,
+ double sppm,
+ double rXs, double rYs, double rZs)
+{
+ double s = sppm * 1e-6;
+ constexpr double radians_per_second = (1.0/(60.0*60.0)) * GPS_PI/180.0;
+ double rX = rXs * radians_per_second; /* radians */
+ double rY = rYs * radians_per_second; /* radians */
+ double rZ = rZs * radians_per_second; /* radians */
+
+ // rotation/scaling matrix = [1+s -rz ry;rz 1+s -rx;-ry rx 1+s];
+ *Dx = tX + (1 + s)*Sx + -rZ*Sy + rY*Sz;
+ *Dy = tY + rZ*Sx + (1 + s)*Sy + -rX*Sz;
+ *Dz = tZ + -rY*Sx + rX*Sy + (1 + s)*Sz;
+}
+
+/* @func GPS_Math_Inverse_Helmert ***********************************
+**
+** Transform from source to coverted datum.
+** translation, scale and rotation parameters are the Helmert parameters
+** to traslate from the converted to source datum!
+**
+** @param [r] Sx [double] source X
+** @param [r] Sy [double] source Y
+** @param [r] Sz [double] source Z
+** @param [w] Dx [double*] converted X
+** @param [w] Dy [double*] converted Y
+** @param [w] Dz [double*] converted Z
+** @param [r] tX [double] translation along x axis (meters)
+** @param [r] tY [double] translation along y axis (meters)
+** @param [r] tZ [double] translation along z axis (meters)
+** @param [r] sppm [double] scale factor - 1 (ppm)
+** @param [r] rXs [double] rotation about x axis (seconds)
+** @param [r] rYs [double] rotation about y axis (seconds)
+** @param [r] rZs [double] rotation about z axis (seconds)
+**
+** @return [void]
+************************************************************************/
+void GPS_Math_Inverse_Helmert(double Sx, double Sy, double Sz,
+ double* Dx, double* Dy, double* Dz,
+ double tX, double tY, double tZ,
+ double sppm,
+ double rXs, double rYs, double rZs)
+{
+ double s = sppm * 1e-6;
+ constexpr double radians_per_second = (1.0/(60.0*60.0)) * GPS_PI/180.0;
+ double rX = rXs * radians_per_second; /* radians */
+ double rY = rYs * radians_per_second; /* radians */
+ double rZ = rZs * radians_per_second; /* radians */
+
+ // forward rotation/scaling matrix is [1+s -rz ry;rz 1+s -rx;-ry rx 1+s]
+ // compute inverse of forward Helmert rotation/scaling matrix
+ // https://www.wolframalpha.com/input?i2d=true&i=Power%5B%7B%7B1%2Bs%2C-Subscript%5Br%2Cz%5D%2CSubscript%5Br%2Cy%5D%7D%2C%7BSubscript%5Br%2Cz%5D%2C1%2Bs%2C-Subscript%5Br%2Cx%5D%7D%2C%7B-Subscript%5Br%2Cy%5D%2CSubscript%5Br%2Cx%5D%2C1%2Bs%7D%7D%2C-1%5D
+ double gain = 1/((s + 1) * (rX*rX + rY*rY + rZ*rZ + s*s + 2*s + 1));
+ double r11 = rX*rX + s*s + 2*s + 1;
+ double r12 = s*rZ + rX*rY + rZ;
+ double r13 = -s*rY + rX*rZ - rY;
+ double r21 = -s*rZ + rX*rY - rZ;
+ double r22 = rY*rY + s*s + 2*s + 1;
+ double r23 = s*rX + rX + rY*rZ;
+ double r31 = s*rY + rX*rZ + rY;
+ double r32 = -s*rX - rX + rY*rZ;
+ double r33 = rZ*rZ + s*s + 2*s + 1;
+
+ *Dx = gain * ((r11 * (Sx-tX)) + (r12 * (Sy-tY)) + (r13 * (Sz-tZ)));
+ *Dy = gain * ((r21 * (Sx-tX)) + (r22 * (Sy-tY)) + (r23 * (Sz-tZ)));
+ *Dz = gain * ((r31 * (Sx-tX)) + (r32 * (Sy-tY)) + (r33 * (Sz-tZ)));
+}
/* @func GPS_Math_Known_Datum_To_WGS84_M **********************************
**
}
+/* @func GPS_Math_WGS84_To_UKOSMap_H ***********************************
+**
+** Convert WGS84 lat/lon to Ordnance survey map code and easting and
+** northing. Uses Helmert.
+**
+** @param [r] lat [double] WGS84 latitude (deg)
+** @param [r] lon [double] WGS84 longitude (deg)
+** @param [w] mE [double *] map easting (metres)
+** @param [w] mN [double *] map northing (metres)
+** @param [w] map [char *] map two letter code
+**
+** @return [int32] success
+************************************************************************/
+int32_t GPS_Math_WGS84_To_UKOSMap_H(double lat, double lon, double* mE,
+ double* mN, char* map)
+{
+ double x;
+ double y;
+ double z;
+ double ax;
+ double ay;
+ double az;
+ double alat;
+ double alon;
+ double ah;
+ double aE;
+ double aN;
+
+
+ GPS_Math_WGS84LatLonH_To_XYZ(lat, lon, 0,
+ &x, &y, &z);
+
+ /* Helmert transformation
+ * https://www.ordnancesurvey.co.uk/documents/resources/guide-coordinate-systems-great-britain.pdf
+ * 6.2 Helmert datum transformations
+ * 6.6 Approximate WGS84 to OSGB36/ODN transformation
+ */
+ // WGS84 -> OSGB36
+ constexpr double tX = -446.448; /* meters */
+ constexpr double tY = +125.157; /* meters */
+ constexpr double tZ = -542.060; /* meters */
+ constexpr double s = +20.4894; /* ppm */
+ constexpr double rX = -0.1502; /* seconds */
+ constexpr double rY = -0.2470; /* seconds */
+ constexpr double rZ = -0.8421; /* seconds */
+
+ GPS_Math_Helmert(x, y, z,
+ &ax, &ay, &az,
+ tX, tY, tZ,
+ s,
+ rX, rY, rZ);
+
+ GPS_Math_XYZ_To_Airy1830LatLonH(&alat, &alon, &ah,
+ ax, ay, az);
+ GPS_Math_Airy1830LatLonToNGEN(alat, alon,
+ &aE, &aN);
+ if (!GPS_Math_EN_To_UKOSNG_Map(aE,aN,mE,mN,map)) {
+ return 0;
+ }
+
+ return 1;
+}
+
+
+/* @func GPS_Math_UKOSMap_To_WGS84_H ***********************************
+**
+** Transform UK Ordnance survey map position to WGS84 lat/lon
+** Uses Helmert transformation
+**
+** @param [r] map [const char *] map two letter code
+** @param [r] mE [double] map easting (metres)
+** @param [r] mN [double] map northing (metres)
+** @param [w] lat [double *] WGS84 latitude (deg)
+** @param [w] lon [double *] WGS84 longitude (deg)
+**
+** @return [int32] success
+************************************************************************/
+int32_t GPS_Math_UKOSMap_To_WGS84_H(const char* map, double mE, double mN,
+ double* lat, double* lon)
+{
+ double E;
+ double N;
+ double alat;
+ double alon;
+ double ht;
+ double ax;
+ double ay;
+ double az;
+ double x;
+ double y;
+ double z;
+
+ if (!GPS_Math_UKOSNG_Map_To_EN(map,mE,mN,&E,&N)) {
+ return 0;
+ }
+
+ GPS_Math_NGENToAiry1830LatLon(E,N,&alat,&alon);
+ GPS_Math_Airy1830LatLonH_To_XYZ(alat, alon, 0.0, &ax, &ay, &az);
+
+ /* Helmert transformation
+ * https://www.ordnancesurvey.co.uk/documents/resources/guide-coordinate-systems-great-britain.pdf
+ * 6.2 Helmert datum transformations
+ * 6.6 Approximate WGS84 to OSGB36/ODN transformation
+ */
+ if constexpr(use_exact_helmert_inverse) {
+ // Actually invert the Helmert transform from WGS84 to OSGB36.
+ // WGS84 -> OSGB36
+ constexpr double tX = -446.448; /* meters */
+ constexpr double tY = +125.157; /* meters */
+ constexpr double tZ = -542.060; /* meters */
+ constexpr double s = +20.4894; /* ppm */
+ constexpr double rX = -0.1502; /* seconds */
+ constexpr double rY = -0.2470; /* seconds */
+ constexpr double rZ = -0.8421; /* seconds */
+
+ GPS_Math_Inverse_Helmert(ax, ay, az,
+ &x, &y, &z,
+ tX, tY, tZ,
+ s,
+ rX, rY, rZ);
+ } else {
+ // Approximate the transform from OSGB36 to WGS84 by using the standard
+ // helmert transform with parameters that all have the opposite signs of
+ // those used to transform from WGS84 to OSGB36.
+ constexpr double tX = +446.448; /* meters */
+ constexpr double tY = -125.157; /* meters */
+ constexpr double tZ = +542.060; /* meters */
+ constexpr double s = -20.4894; /* ppm */
+ constexpr double rX = +0.1502; /* seconds */
+ constexpr double rY = +0.2470; /* seconds */
+ constexpr double rZ = +0.8421; /* seconds */
+
+ GPS_Math_Helmert(ax, ay, az,
+ &x, &y, &z,
+ tX, tY, tZ,
+ s,
+ rX, rY, rZ);
+ }
+
+ GPS_Math_XYZ_To_WGS84LatLonH(lat, lon, &ht, x, y, z);
+
+ return 1;
+}
+
/* @func GPS_Math_WGS84_To_UKOSMap_C ***********************************
**
double Sif, double* Dphi, double* Dlam,
double* DH, double Da, double Dif, double dx,
double dy, double dz);
+void GPS_Math_Helmert(double Sx, double Sy, double Sz,
+ double* Dx, double* Dy, double* Dz,
+ double tX, double tY, double tZ,
+ double sppm,
+ double rXs, double rYs, double rZs);
+void GPS_Math_Inverse_Helmert(double Sx, double Sy, double Sz,
+ double* Dx, double* Dy, double* Dz,
+ double tX, double tY, double tZ,
+ double sppm,
+ double rXs, double rYs, double rZs);
void GPS_Math_Known_Datum_To_WGS84_M(double Sphi, double Slam, double SH,
double* Dphi, double* Dlam, double* DH,
int32_t n);
double* mN, char* map);
int32_t GPS_Math_UKOSMap_To_WGS84_M(const char* map, double mE, double mN,
double* lat, double* lon);
+int32_t GPS_Math_WGS84_To_UKOSMap_H(double lat, double lon, double* mE,
+ double* mN, char* map);
+int32_t GPS_Math_UKOSMap_To_WGS84_H(const char* map, double mE, double mN,
+ double* lat, double* lon);
int32_t GPS_Math_WGS84_To_UKOSMap_C(double lat, double lon, double* mE,
double* mN, char* map);
int32_t GPS_Math_UKOSMap_To_WGS84_C(const char* map, double mE, double mN,
&result);
valid = (ct == 3);
if (valid) {
- if (! GPS_Math_UKOSMap_To_WGS84_M(map, lx, ly, &lat, &lon))
+ if (! GPS_Math_UKOSMap_To_WGS84_H(map, lx, ly, &lat, &lon))
fatal("%s: Unable to convert BNG coordinates (%s)!\n",
module, str);
}
--- /dev/null
+No,BNG-Zone,BNG-East,BNG-North,Name,Description,Date,Time\r
+1,ST,96552, 3097,"Badbury Rings Trig Point","Test Point 1",2024/08/05,00:00:00\r
+2,ST,91437, 1894,"Spettisbury Rings Trig Point","Test Point 2",2024/08/05,00:00:00\r
--- /dev/null
+<?xml version="1.0" encoding="UTF-8"?>
+<gpx version="1.0" creator="GPSBabel - https://www.gpsbabel.org" xmlns="http://www.topografix.com/GPX/1/0">
+ <time>1970-01-01T00:00:00Z</time>
+ <bounds minlat="50.816527451" minlon="-2.122918364" maxlat="50.827399503" maxlon="-2.050323375"/>
+ <wpt lat="50.827399503" lon="-2.050323375">
+ <time>2024-08-05T00:00:00Z</time>
+ <name>Badbury Rings Trig Point</name>
+ <cmt>Test Point 1</cmt>
+ <desc>Test Point 1</desc>
+ </wpt>
+ <wpt lat="50.816527451" lon="-2.122918364">
+ <time>2024-08-05T00:00:00Z</time>
+ <name>Spettisbury Rings Trig Point</name>
+ <cmt>Test Point 2</cmt>
+ <desc>Test Point 2</desc>
+ </wpt>
+</gpx>
<?xml version="1.0" encoding="UTF-8"?>
<gpx version="1.0" creator="GPSBabel - https://www.gpsbabel.org" xmlns="http://www.topografix.com/GPX/1/0">
<time>1970-01-01T00:00:00Z</time>
- <bounds minlat="50.207517638" minlon="-7.489235772" maxlat="60.150151215" maxlon="1.296391371"/>
- <wpt lat="58.007809302" lon="-6.751777637">
+ <bounds minlat="50.207587993" minlon="-7.489131626" maxlat="60.150023344" maxlon="1.296336192"/>
+ <wpt lat="58.007715410" lon="-6.751692411">
<time>2008-08-23T20:56:12Z</time>
<name>AlineLodge</name>
<cmt>Aline Lodge</cmt>
<desc>Aline Lodge</desc>
<sym>City (Small)</sym>
</wpt>
- <wpt lat="53.140870076" lon="-4.273683778">
+ <wpt lat="53.140881274" lon="-4.273629535">
<time>2008-08-23T20:55:54Z</time>
<name>Caernarfon</name>
<cmt>Caernarfon</cmt>
<desc>Caernarfon</desc>
<sym>City (Small)</sym>
</wpt>
- <wpt lat="50.207517638" lon="-5.295409027">
+ <wpt lat="50.207587993" lon="-5.295322139">
<time>2008-08-23T20:55:39Z</time>
<name>Camborne</name>
<cmt>Camborne</cmt>
<desc>Camborne</desc>
<sym>City (Small)</sym>
</wpt>
- <wpt lat="50.712891291" lon="-3.537602671">
+ <wpt lat="50.712954134" lon="-3.537551739">
<time>2008-08-23T20:55:35Z</time>
<name>Exeter</name>
<cmt>Exeter</cmt>
<desc>Exeter</desc>
<sym>City (Medium)</sym>
</wpt>
- <wpt lat="56.634525358" lon="-2.889408436">
+ <wpt lat="56.634466564" lon="-2.889403757">
<time>2008-08-23T20:55:22Z</time>
<name>Forfar</name>
<cmt>Forfar</cmt>
<desc>Forfar</desc>
<sym>City (Small)</sym>
</wpt>
- <wpt lat="55.426025761" lon="-2.790531785">
+ <wpt lat="55.425992284" lon="-2.790521406">
<time>2008-08-23T20:55:25Z</time>
<name>Hawick</name>
<cmt>Hawick</cmt>
<desc>Hawick</desc>
<sym>City (Small)</sym>
</wpt>
- <wpt lat="57.621914847" lon="-7.489235772">
+ <wpt lat="57.621828064" lon="-7.489131626">
<time>2008-08-23T20:56:17Z</time>
<name>Hosta</name>
<cmt>Hosta</cmt>
<desc>Hosta</desc>
<sym>City (Small)</sym>
</wpt>
- <wpt lat="57.007373614" lon="-6.271139683">
+ <wpt lat="57.007301216" lon="-6.271060826">
<time>2008-08-23T20:56:24Z</time>
<name>IsleofRhum</name>
<cmt>Isle of Rhum</cmt>
<desc>Isle of Rhum</desc>
<sym>City (Small)</sym>
</wpt>
- <wpt lat="60.150151215" lon="-1.142586363">
+ <wpt lat="60.150023344" lon="-1.142651285">
<time>2008-08-23T20:59:39Z</time>
<name>Lerwick</name>
<cmt>Lerwick</cmt>
<desc>Lerwick</desc>
<sym>City (Small)</sym>
</wpt>
- <wpt lat="58.684841865" lon="-3.116680165">
+ <wpt lat="58.684739998" lon="-3.116684802">
<time>2008-08-23T20:57:48Z</time>
<name>Nethertown</name>
<cmt>Nethertown</cmt>
<desc>Nethertown</desc>
<sym>City (Small)</sym>
</wpt>
- <wpt lat="52.635495464" lon="1.296391371">
+ <wpt lat="52.635529551" lon="1.296336192">
<time>2008-08-23T20:55:44Z</time>
<name>Norwich</name>
<cmt>Norwich</cmt>
<desc>Norwich</desc>
<sym>City (Medium)</sym>
</wpt>
- <wpt lat="50.741731984" lon="-0.782771513">
+ <wpt lat="50.741799936" lon="-0.782773783">
<time>2008-08-23T20:57:02Z</time>
<name>Selsey</name>
<cmt>Selsey</cmt>
<desc>Selsey</desc>
<sym>City (Small)</sym>
</wpt>
- <wpt lat="53.393555571" lon="-1.494139002">
+ <wpt lat="53.393567055" lon="-1.494143092">
<time>2008-08-23T20:55:01Z</time>
<name>Sheffield</name>
<cmt>Sheffield</cmt>
gpsbabel -i unicsv -f ${REFERENCE}/unidelim.csv -o gpx -F ${TMPDIR}/unidelim.gpx
compare ${REFERENCE}/unidelim.gpx ${TMPDIR}/unidelim.gpx
+gpsbabel -i unicsv,utc -f ${REFERENCE}/bngtest.csv -o gpx -F ${TMPDIR}/bngtest~csv.gpx
+compare ${REFERENCE}/bngtest.gpx ${TMPDIR}/bngtest~csv.gpx
+gpsbabel -i gpx -f ${REFERENCE}/bngtest.gpx -o unicsv,grid=bng,utc -F ${TMPDIR}/bngtest~gpx.csv
+compare ${REFERENCE}/bngtest.csv ${TMPDIR}/bngtest~gpx.csv
+
#include "garmin_fs.h" // for garmin_fs_t
#include "garmin_tables.h" // for gt_lookup_datum_index, gt_get_mps_grid_longname, gt_lookup_grid_type
#include "geocache.h" // for Geocache, Geocache::status_t, Geoc...
-#include "jeeps/gpsmath.h" // for GPS_Math_UKOSMap_To_WGS84_M, GPS_Math_EN_To_UKOSNG_Map, GPS_Math_Known_Datum_To_UTM_EN, GPS_Math_Known_Datum_To_WGS84_M, GPS_Math_Swiss_EN_To_WGS84, GPS_Math_UTM_EN_To_Known_Datum, GPS_Math_WGS84_To_Known_Datum_M, GPS_Math_WGS84_To_Swiss_EN, GPS_Math_WGS...
+#include "jeeps/gpsmath.h" // for GPS_Math_UKOSMap_To_WGS84_H, GPS_Math_EN_To_UKOSNG_Map, GPS_Math_Known_Datum_To_UTM_EN, GPS_Math_Known_Datum_To_WGS84_M, GPS_Math_Swiss_EN_To_WGS84, GPS_Math_UTM_EN_To_Known_Datum, GPS_Math_WGS84_To_Known_Datum_M, GPS_Math_WGS84_To_Swiss_EN, GPS_Math_WGS...
#include "session.h" // for session_t
#include "src/core/datetime.h" // for DateTime
#include "src/core/logging.h" // for Warning, Fatal
fatal(MYNAME ": Unable to convert BNG coordinates (%.f %.f)!\n",
bng_easting, bng_northing);
}
- if (! GPS_Math_UKOSMap_To_WGS84_M(
+ if (! GPS_Math_UKOSMap_To_WGS84_H(
bngz, bnge, bngn,
&wpt->latitude, &wpt->longitude))
fatal(MYNAME ": Unable to convert BNG coordinates (%s %.f %.f)!\n",
bngz, bnge, bngn);
} else { // traditional zone easting northing
- if (! GPS_Math_UKOSMap_To_WGS84_M(
+ if (! GPS_Math_UKOSMap_To_WGS84_H(
CSTR(bng_zone), bng_easting, bng_northing,
&wpt->latitude, &wpt->longitude))
fatal(MYNAME ": Unable to convert BNG coordinates (%s %.f %.f)!\n",
double north;
double east;
- if (! GPS_Math_WGS84_To_UKOSMap_M(wpt->latitude, wpt->longitude, &east, &north, map)) {
+ if (! GPS_Math_WGS84_To_UKOSMap_H(wpt->latitude, wpt->longitude, &east, &north, map)) {
unicsv_fatal_outside(wpt);
}
auto fieldWidth = fout->fieldWidth();
#include "garmin_fs.h" // for garmin_fs_t
#include "geocache.h" // for Geocache, Geocache::status_t, Geoc...
#include "grtcirc.h" // for RAD, gcdist, radtometers
-#include "jeeps/gpsmath.h" // for GPS_Math_WGS84_To_UTM_EN, GPS_Lookup_Datum_Index, GPS_Math_Known_Datum_To_WGS84_M, GPS_Math_UTM_EN_To_Known_Datum, GPS_Math_WGS84_To_Known_Datum_M, GPS_Math_WGS84_To_UKOSMap_M
+#include "jeeps/gpsmath.h" // for GPS_Math_WGS84_To_UTM_EN, GPS_Lookup_Datum_Index, GPS_Math_Known_Datum_To_WGS84_M, GPS_Math_UTM_EN_To_Known_Datum, GPS_Math_WGS84_To_Known_Datum_M, GPS_Math_WGS84_To_UKOSMap_H
#include "jeeps/gpsport.h" // for int32
#include "session.h" // for session_t
#include "src/core/datetime.h" // for DateTime
char map[3];
double north;
double east;
- if (! GPS_Math_WGS84_To_UKOSMap_M(wpt->latitude, wpt->longitude, &east, &north, map))
+ if (! GPS_Math_WGS84_To_UKOSMap_H(wpt->latitude, wpt->longitude, &east, &north, map))
fatal(MYNAME ": Position (%.5f/%.5f) outside of BNG.\n",
wpt->latitude, wpt->longitude);
buff = QString::asprintf(fmp.printfc.constData(), map, qRound(east), qRound(north));
projects / gpsbabel.git / commitdiff